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Missile  autupilol  design  which  iiicurpiirales  a digital  eumputer  iir  inieruprueessorls)  to  implement  the 
stability  compensation  equations  requires  considerations  over  and  above  those  made  for  analog  computer 
autopilots.  Some  of  these  cuiisideratioiis  are:  the  selection  of  sampling  frequency,  the  A D converter 
wordlength,  the  phase  lag  contributed  by  the  zero  order  data  hold  (i.e.,  the  D A converter),  the  transport 
lag  in  the  digital  computer,  and  the  effects  of  sampling  high  frequency  sensor  and  structural  (bending)  noise 
(foldover  or  “aliasing").  Ttie  purpose  of  this  paper  is  to  discuss  some  considerations  in  the  selection  of  a 
sample  rate  for  the  digital  autopilot. 


In  the  following  discussion  a distinction  is  made  between  guidance  and  control.  Guidance  is  assumed 
to  be  that  aspect  of  the  missile  which  determines  the  required  flight  path  of  the  missile  which  causes  target 
intercept.  Control  is  assumed  to  be  that  portion  which  provides  satisfactory  missile  behavior  about  that 
flight  path. 


ANALOG  VKRSUS  DIGITAL  CONTROL  COMPENSATION  (AUTOPILOT) 


The  analog  control  computer  processes  the  error  sensor  outputs  by  comparing  the  measured  response 
to  the  required  or  comtiunded  resporse,  and  alters  the  error  signal  to  provide  correction  commands  to  give 
desired  response  without  oscillation. 


The  digital  control  computer  also  makes  comparisons  of  desired  versus  measured  response  to 
determine  the  correct  response  command.  However,  the  error  determination  is  corrupted  by  noninfinite 
resolution  A.'D  converters.  The  real  world  AID  converter  not  only  gives  quantization  noise  but  often 
generates  bit  noise  which  further  corrupts  the  error  signals.  At  the  opposite  end  of  the  computer,  an 
additional  problem  is  encountered.  The  digital  computer  can  only  determine  the  error  signal  once  every 
sample;  therefore,  the  output  is  not  smooth  but  is  constant  (with  a zero  order  hold)  between  samples.  This 
has  the  effect  of  adding  a delay  to  the  system,  and  a delay  is  almost  always  destabilizing.  Thirdly,  the 
computer  requires  time  to  process  the  input  signals  in  order  to  generate  the  output  command.  That  time 
required  is  a second  transport  lag  that,  in  general,  is  destabilizing. 


If  a satisfactory  analog  control  autopilot  has  been  formulated  and  tested,  and  then  a digital  autopilot 
is  proposed,  a common  approach  is  to  emulate  the  analog  autopilot.  This  is  especially  common  in  situations 
which  the  plant  (missile  plus  actuator,  etc.)  is  not  totally  known  because  of  lost  aerodynamic  data,  actuator 
data,  or  whatever.  This  is  the  cause  of  much  concern  by  the  digital  autopilot  designer  because  in  this  case  he 
does  not  really  know  just  how  much  gain  and  phase  margins  existed  in  the  first  place.  Thus  he  is  forced  into 
a situation  where  he  tries  to  maintain  the  same  gain  and  phase  shift  through  the  digital  computer  as  was 
specified  for  the  analog  control  computer.  If  he  succeeds  perfectly  in  gain  and  phase  relationships,  he  still  is 
burdened  with  the  noise  introduced  by  the  sampler  (A-D  converter).  In  many  missile  systems,  additional 
noise  injected  into  the  control  loop  has  extremely  detrimental  effects  on  the  actuation  system,  usually  in 
the  form  of  much  greater  hydraulic  or  electric  power  requirements.  Thus  the  control  engineer  is  faced  with 
the  requirement  to  emulate  an  analog  system  (which  cannot  be  perfectly  accomplished)  without  utilizing 
the  advantages  of  flexibility  and  versatility  offered  by  the  digital  autopilot.  Flexibility  allows  the 
compensation  to  be  changed  almost  constantly  during  flight  to  provide  optimum  stability  margins  rather 
than  just  one  or  two  “bands"  which  are  a compromise  between  the  operating  extremes  expected. 
Versatility  enables  the  system  designer  to  incorporate  options  not  readily  available  to  the  analog  designer, 
such  as  a wide  variety  of  open  loop  maneuvers,  and  to  change  the  guidance  law  for  different  portions  of 
flight. 
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DIGITAL  CONTROL  SYSTEM  DESIGN 
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In  order  to  utilize  the  analog  control  engineer  to  the  maximum  extent  in  the  design  of  digital  control 
systems,  one  technique  often  used  is  to  develop  an  analog  control  compensator  and  then  use  the  “Kader 
and  Gold”  method  (Kef.  1)  of  converting  the  compensation  to  a digital  compensator.  This  method 
preserves  the  gain  and  phase  characteristics  of  the  analog  compensator  to  a frequency  up  to  approximately 
one-third  the  sample  frequency.  It  consists  of  prewtrping  the  break  frequencies  and  then  making  the 
following  substitution  for  s 


s - 


z-1 


The  prewarping  consists  of  altering  the  break  frequency  (a)  by  the  following  equation  , eT  where 

A " mH  2 

T - sample  period.  This  design  method  works  reasonably  well.  The  control  engineer  frequently  selects  a 
sample  frequency  ten  times  the  frequency  of  interest  in  order  to  maintain  the  gain  and  phase  margins 
determined  in  his  analog  design.  This  method  however,  does  not  account  for  the  transport  lag  of  the  digital 
computer  nor  the  zero  order  hold  phase  lag  unless  approximations  of  these  terms  had  been  included  in  the 
analog  system  design.  By  not  including  these  terms,  the  phase  margins  suffer  by  the  following  amounts.  The 
zero  order  hold  phase  lag  is 


180 


0. 


LAG(ZOH) 


(degrees) 


where  f^  is  the  critical  frequency  (usually  rigid  body  crossover  frequency)  and  is  the  sample  frequency. 


Thuiifthcdesignbadbecnielected  where  there  was  a 30”  phase  margin  at  the  Odh  crossover  frequency, 
the  digital  version  would  have  only  12''  phase  margin  assuming  the  sample  rate  10  times  crossover 
frequency  had  been  selected.  This  phase  margin  is  unacceptably  small,  si  it  means  greater  than  a 14  dB 
resonant  peak  in  the  closed  loop  frequency  response. 


In  addition,  the  computational  delay  in  the  digital  airborne  computer  also  provides  a phase  lag  of  the 
following  amount 


where  t|j  - delay  in  seconds 


^LAG(DELAY) 


- 360  f^  X t|j  (degrees) 


In  practice,  the  digital  airborne  computer  is  sufficiently  loaded  as  to  require  between  50%  and  80%  of 
its  sample  period  to  sample  the  input  error  senson  and  compute  the  required  output  commands.  Thus  the 
delay  can  be  found  to  be  within  the  following  range 


360* 


0.5 


0.8 


•fr  < 0 J V < 360*  - . 

® delay  (degrees)  f 
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If  the  sample  frequency  had  been  selected  to  be  ten  times  crossover  frequency,  then  the  delay  phase  lag  has 
the  range 


18°  < « < 28.8°. 
LAG(DELAY) 
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Thi'  total  phase  lag  for  this  typical  case,  with  both  zero  order  hold  and  computational  delay,  is  between  kc” 
and  46.8  which  can  give  a negative  phase  margin  m the  system  if  these  effects  had  not  been  included  in  the 
original  analog  design.  To  have  included  them  in  the  original  analog  design  wiiuld  have  required  the  use  of 
I’ade  approximations  ( Kef.  2)  which  increase  the  order  of  the  system  by  the  number  of  terms  chosen  in  the 
approximation. 

A mure  direct  approach  to  digital  autopilot  design  is  through  the  use  of  / transforms.  This  technique 
involves  no  approximations  of  transport  lags  or  zero  order  holds.  Then'  are  several  / transform  computer 
programs  iKef.  3|  which  utilize  the  theory  of  residues  in  order  to  avoid  series  truncations  to  find 
/transforms. 

These  programs  accx'pt  the  t.aplace  transform  definition  of  the  plant  and  computational  delay  and 
give  a ratio  of  polynominals  in  /-domain,  for  which  a suitable  compimsation  is  to  be  determined.  This 
technique  does  not  tell  the  designer  what  sample  rate  to  sided,  but  oy  choosing  different  sample  rates  and 
compensating  the  resultant  equations,  the  designer  can  determine  which  sample  rate  best  allows  him  to 
meet  stability  criteria  with  the  minimum  penalty  (compensator  order  and  amount  of  “lead"  required).  It 
has  been  shown  that  with  missiles  having  lightly  damped  body  bending  roles,  a high  sample  rate  is  most 
lieneficial  (Ref.  4).  This  is  true  because  the  phase  lags  of  the  computational  delay  and  zero  order  hold  are 
still  present  and  must  be  accounted  for  by  sufricient  phase  lead  in  the  compensation.  These  items  have  not 
disappeared  but  have  become  more  directly  asses.sable  by  the  / transform  technique. 

In  general,  analog  control  engineers  are  reluctant  to  use  the  /-transform  fur  the  following  reasons, 
first,  a /-transform  program  is  nut  usually  immediately  available,  or  if  It  is  available,  there  is  no  confidence 
in  it  because  of  the  lack  of  experience.  Secondly  it  requires  a different  criteria  fur  stability  in  that  the  roots 
of  the  open  loop  transfer  function  must  have  magnitude  less  than  unity  (quadratic  routs  must  be  in  polar 
form)  and  stability  margins  in  the  /-domain  are  somewhat  more  difficult  to  assess,  i.e.,  lines  of  constant 
damping  in  the  unit  circle  are  distorted. 


In  order  to  get  around  these  objectives  and  to  let  the  analog  control  engineer  use  continuous  system 
analysis  and  design  tools,  a transformation  has  been  made  available  which  maps  the  /-domain  unit  circle 
into  a left  half  plane  of  a new  variable  W.  The  change  of  variables  is 


and  is  a by-product  of  some  /-trans'orm  progranxs  (Ref.  3).  By  the  use  of  this  transformation,  the 
compensation  designer  can  treat  the  op?n  loop  transfer  function  in  W-domain  as  if  it  were  in  the  S-domain, 
and  all  programs  for  root  locus,  Bode  p ots,  Nichols  charts,  Nyquist  plots,  etc.  for  continuous  systems  are 
useable.  Again,  this  does  not  reduce  the  problems  associated  with  digital  control  autopilots  but  merely 
makes  them  more  addressable, 

A number  of  digital  autopilot  anilysts  h,ive  attempted  to  utilize  sumulation  techniques  in 
transforming  the  analog  design  to  a digital  design.  Others  have  tried  matching  the  impulse  response  of  the 
analog  compensator  by  a difference  equation  having  the  same  impulse  response.  It  has  been  noted  that 
those  two  methods  fail  to  account  for  the  computer  delay  and  output  hold,  but  also,  they  do  not  match  the 
zero  (numerator)  frequency  response  of  the  desired  compensator.  (See  appendix  A).  In  general,  simulation 
techniques  should  be  avoided  in  sampled-data  autopilot  design. 
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BOIJY  BENDING 


Su  far  only  ri|jid  body  crossover  frequencies  and  the  phase  lag  of  the  digital  system  at  these 
frequencies  have  been  discussed.  The  body  bending  of  missiles  cannot  be  dismissed  because  quite 
often  the  actuation  (controller)  technique  is  sensitive  to  signals  of  these  frequencies.  Blowdown 
type  actuators  will  require  rnoiv  res»Tve  if  the  actuator  encounters  the  body  bending  frequency. 
Hydraulic  actuators  may  require  greater  pump  capacity  or  accumulator  reserve  if  these  frequencies 
exist  for  any  length  of  time.  In  a relatively  large  percent  of  the  cases,  the  amplitude  of  the 
missile  body  bending  (mode  1)  is  sufficiently  great  to  cause  the  loop  gain  to  be  greater  than  0 
dB,  therefore  potentially  unstable.  The  designer  cannot  always  reduce  this  amplitude  sufficiently, 
either  because  of  its  proximity  to  rigid  body  crossover  frequency  or  because  its  range  is  too  broad 
for  a notch  filter.  Therefore,  this  mode  is  "phase  stabilized’’  by  making  the  phase  at  that 
frequency  range  sufficiently  far  from  180°  that  instability  cannot  occur.  (See  Figure  1).  In  analog 
systems  a "lead”  type  compensator  provides  just  enough  phase  lead  to  give  the  required  pha.se 
margin  (there  is  no  gain  margin).  .Should  an  attempt  to  "digitize"  this  analog  design  be  made 
without  accounting  fur  the  computational  and  hold  lags,  the  phase  margin  would  be  lost  with 
tragic  results.  For  example.  If  mode  1 body  bending  had  a frequency  of  & times  rigid  body 
crossover  frequency,  then  the  phase  lag  due  to  the  delay  and  hold  would  be 

<^delay 

where  sample  frequency  is  10  times  rigid  body  crossover  frequency.  For  that  reason,  many  analog  control 
engineers  who  are  attempting  to  phase-stabilize  mode  1 body  bending  will  select  a sample  rate  of  10-15 
times  the  bending  frequency,  thereby  assuring  that  the  delays  of  20°  to  30°  can  be  relatively  easily 
compensated  for  in  in  the  compensation  design.  ’The  same  thing  must  be  done  in  the  direct  Z or  \tmomain, 
but  in  a more  direct  fashion.  In  an  analytical  sense,  the  designer  can  provide  an  extremely  large  amount  of 
lead,  including  multiple  pure  zeroes,  but  practical  experience  shows  that  this  kind  of  design  would  gicatly 
accentuate  any  sensor  and  A-0  converter  noise  would  quickly  deplete  the  actuation  system.  By  selecting 
high  sample  frequency  to  minimize  the  amount  of  lead  compensation  required  to  account  for  the  delay,  the 
autopilot  designer  suddenly  realizes  that  he  has  chosen  a sample  which  operates  at  10-15  times  mode  1 
body  bending  frequency,  which  is  perhaps  50-75  times  rigid  body  crossover  frequency  - and  this  may  be 
100  to  150  times  the  missile  rigid  body  airframe  natural  frequency. 

These  high  sample  frequencies  introduce  additional  penalties.  It  has  been  shown  (Ref.  1,  p.  102)  that 
the  high  sample  rates  generate  new  problems  such  as  greater  required  minimum  A/D  converter  wordlength 
(difricult  to  obtain  at  high  sample  rates)  and  internal  computer  wordlength  or  use  of  double  precision 
arithmetic  which  increases  computational  delay. 

Signal  foldover,  also  called  aliasing,  has  been  given  considerable  attention  by  authors  of  digital 
autopilot  theory  papers  (Ref.  5).  In  light  of  the  high  sample  rate  selection  for  phase  considerations,  signal 
foldover  of  rigid  body  and  mode  1 bending  is  not  a problem.  Sensor  noise  near  the  sample  frequency  will 
probably  be  of  low  amplitude  and  folded  signals  should  then  be  of  little  consequence.  However,  for  safety 
margin,  the  compensation  should  include  a second  order  lag  filter  with  a break  frequency  just  high  enough 
to  have  little  effect  on  the  phase  lag  at  the  frequency  of  Interest,  which  is  rigid  body  crossover  of  mode  1 
bending.  Communication  type  niters  such  as  Butterworth,  Chebycheff,  etc.,  cannot  be  used  because  of  the 
large  phase  lags  associated  with  those  niters  at  low  frequency  compared  to  their  stopbands. 
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CAIM  MARGIN 


Figure  1.  Nichols  chart  for  sampicii-data  autopiht. 


orn.MAI.  COM  ROLLERS 


Optimal  coiitrul  llu-ory  has  bum  applied  to  digital  autopilots  'Ref.  6)  but,  in  general,  this  technique 
does  not  offer  some  of  the  advantages  which  clas.sieal  theor>'  provides.  First,  optimal  control  theory 
assumes  that  the  system  is  a regulator  and  that  the  designer  is  attempting  to  ivro  all  states  with  some 
miimni/ation  of  control  or  stale  action.  It  alsoassuim-s  infinite  time  to  accomplish  this  result.  The  problem 
is  presented  m this  manner  to  keep  it  mathematically  tractable.  Unfortunately,  the  feedback  matrix 
coefficients  are  strongly  dependent  on  knowledge  of  the  plant  (the  A matrix)  and  feedback  in  the  classical 
sense  to  compensate  fur  uncertainties  in  the  plant  parameters  is  not  provided.  All  system  slates  must  be  fed 
back  III  the  optimal  controller.  Unmeasureable  stales  must  be  synthesized  by  stale  observers  or  estimators, 
which  further  increase  the  complexity  of  the  total  digital  control  autopilot.  Although  optimal  control 
theory  does  not  provide  any  quantum  jump  in  stability  margins  over  classical  control  theory,  it  does 
provide  for  controlling  several  variables,  which  is  an  improvement  over  single-input  single-output  systems. 

.Stochastic  optimal  control  theory,  in  the  form  of  Kalman  filter  theory,  has  provided  results  to 
estimate  plant  stales  from  signals  corrupted  by  noise,  which  has  proven  successful  in  the  SAM-D  autopilot. 
Digital  autopilot  design  should  certainly  consider  the  Kalman  filter,  although  it  falls  more  in  the  domain  of 
guidance  than  control.  It  should  be  considered  then  in  the  decision  of  digital  computer  sizing  and  memory 
as  it  offers  benefits  unavailable  in  most  analog  airborne  computers.  The  guidance  portion  of  the  guidance 
and  control  of  missiles  is  relatively  low  frequency,  on  the  order  of  a tenth  of  the  rigid  body  control 
frequencies.  Con.sequenlly,  the  guidance  equation  is  solved  at  a much  slower  rate  than  control,  and  only  a 
portion  of  the  guidance  equations  are  solved  in  each  control  sample  period.  It  may  take  5 to  10  control 
sample  periods  to  output  a guidance  command.  This  technique  offers  a problem  for  the  control  engineer,  in 
that  the  entire  guidance  and  control  system  is  now  a multiloop,  multirate  system.  In  the  Z-domain 
comjiensation  design  this  problem  is  frequently  not  solvable,  as  block  diagram  reduction  for  nonerror 
sampled  systems  is  nut  possible  ( Ref  7 ).  The  autopilot  designer  must  then  assume  that  loop  interactions  do 
not  occur. 

CONCLUSION 


This  report  has  discussed  some  of  the  considerations  in  the  design  of  digital  autopilots  for  missiles.  A 
major  aspect  of  the  design  is  the  selection  of  the  sampling  rate,  since  this  factor  dictates  A-D  converter 
rates,  wordlength,  and  the  computer  speed  requirement.  Because  the  computational  delay  and  the  zero 
order  hold  impart  phase  lag  inversely  proportional  to  sampling  frequency,  the  designer  mutt  select  a sample 
rate  10  to  lb  times  mode  1 body  bending  frequency  if  it  is  phase  stabilized.  If  bending  amplitude  is 
sufficiently  low,  a sample  rate  selection  of  10  to  15  rigid  body  crossover  frequency  is  sufficient  to  assure 
adequate  rigid  body  phase  margin.  The  digital  computer  provides  the  ability  to  adjust  the  compensation 
throughout  the  flight  to  obtain  maximum  vehicle  response  within  stability  constraints,  while  most  analog 
autopilots  have  a single  configuration  which  is  required  to  give  reasonable  stability  margin  over  a wide  range 
of  flight  conditions.  The  guidance  aspect  of  missile  guidance  and  control  can  most  benefit  from  the  modern 
techniques  of  Kalman  filteri  ,g  and  optimal  response. 
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APPENDIX  A 

IMPULSE  INVARIANT  TRANSFORMATIONS 
FOR  DIGITAL  COMPENSATION 


Thr  respunsr  uf  an  initially  relaxed  system  of  differential  equations  to  a unit  impulse  is 
equal  to  the  inverse  Laplace  transform  of  the  transfer  function.  Thus  the  transfer  function  is 
really  the  transform  of  the  impulse  response. 

If  the  impulse  response  uf  the  system  is  sampled,  the  Z transform  of  the  respon.se  gives  the 
output  at  the  sample  instants.  The  Z-transfurm  is  referred  to  as  an  impulse-invaiiant  transformation 
becau.se  of  that  characteristic.  The  following  diagram  depicts  the  situation. 


R(l)-S(t) 


•(s).C(z) 


where  G(z)  - Z ( G*(s)) 


R(«)  - 1 


If  G(s)  is  an  integrator  — , the  following  output  results: 

I 

C(z)  - - Z ( G*(»)  )• 

Hie  inverse  of  this  transform  gives  the  time  function  of  the  output.  In  this  case, 

C»(t)  - 1, 

which  means  an  output  of  unit  amplitude  at  each  sample  instant. 

Control  Engineers  seek  compensation  equations  which  are  transfer  functions  whose  impulse 
response  causes  the  system  to  meet  error  and  stability  requirements.  Some  engineen  have 
attempted  to  “digitize"  an  analog  compensator  by  implementing  the  Z-transform  by  a set  of 
difference  equations.  They  feel  that  since  the  Z-transform  is  an  impulse  invariant  transfoimation, 
this  technique  should  work.  However,  that  method  is  incorrect  for  the  following  reason.  The 
replacement  of  the  analog  transfer  function  by  a difference  equation  set  is  shown  in  the  following 
diagram. 


The  situation  is  dramatically  different  fi’om  the  preceeding  diagram  by  the  addition  of  the  samples 
in  front  of  the  compensator.  As  is  pointed  out  in  many  tests  (1), 


( R(s)G(s)  ) - RG(z)  ,t  R(z)G(z). 


(1)  “Introduction  to  Continuous  and  Digital  Control  Systems.”  Roberts  Saucedo  and  Earle 
Schiring,  MacMillian,  1968,  p.  231. 
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' In  thr  ^i^^t  ri|ure,  the  iiulput  was 

IC(z)  - RG(z); 
in  the  second  figure 

C(z)  - R(z)0(z|. 

A few  trials  with  various  inputs  K(Z)  will  demonstrate  the  inequality.  This  problem  is  not 
^ dependent  on  sampie  rate,  but  is  always  true. 

^1  The  same  sort  of  matching  technique  described  above  has  been  tried  by  engineers  using  the 

' State  Space  methods  of  .Modern  Control  Theory.  In  their  method  the  analog  compensation 

equations  are  rewritten  as  first  order  differential  equations  called  .State  Equations.  Then  they  attempt  to 
I find  a set  of  first  order  difference  equations  whose  outputs  at  the  sample  Instants  It  equal  to  the  analog 

state  equations  excited  by  a unit  impulse  Input  vector.  That  is  equivalent  to  the  flrat  described  classical 
[ approach  and  equally  erroneous  results  will  be  obtained. 

I 

■ j 
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